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Abstract 

We investigate polydispersity effects on tlie average structure factor of colloidal 
suspensions of neutral particles with surface adhesion. A sticky hard sphere model 
alternative to Baxter's one is considered. The choice of factorizable stickiness param- 
eters in the potential allows a simple analytic solution, within the "mean spherical 
approximation", for any number of components and arbitrary stickiness distribu- 
tion. Two particular cases are discussed: i) all particles have different sizes but equal 
stickiness (Model I), and ii) each particle has a stickiness proportional to its size 
(Model II). The interplay between attraction and polydispersity yields a markedly 
different behaviour for the two Models in regimes of strong coupling (i.e. strong ad- 
hesive forces and low temperature) and large polydispersity. These results are then 
exploited to reanalyze experimental scattering data on sterically stabilized silica 
particles. 

Keywords: Colloidal models, sticky hard spheres, polydispersity, structure factors, 
small angle scattering. 



1 Introduction 



In studies on colloidal suspensions of neutral particles with adhesive inter- 
actions [1-3], several attempts have been made to fit small angle scattering 
data by using Baxter's "sticky hard sphere" (SHS) model [4]. In addition to 
a hard sphere (HS) repulsion, the SHS potential contains a surface adhesive 
term, introduced by Baxter as a limit of an attractive square well which be- 
comes infinitly deep and narrow according to a particular prescription [4] . For 
this model the Ornstein-Zernike (OZ) integral equations of the liquid state 
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theory have been solved analytically within the Percus-Yevick (PY) approxi- 
mation, and expressions for the correlation functions are thus available [4-6]. 
Unfortunately, however, for a fluid with p components the SHS-PY analytical 
solution requires the knowledge of a set of density- dependent parameters Xij, 
to be determined through p(p + l)/2 coupled quadratic equations [5,6]. This 
feature makes Baxter's SHS-PY solution for mixtures practically inapplicable 
to colloidal systems with significant polydispersity, i.e. with large p. 

Polydispersity means that the mesoscopic suspended particles of a same chem- 
ical species may differ in size, charge or other properties. Consequently, a 
polydisperse fluid is always a mixture with a large number of components and 
strong asymmetries (in size, charge, etc.), even when all meso-particles belong 
to a unique chemical species. 

In a recent paper (hereafter referred to as Reference I) [7] we showed that 
the aforesaid shortcomings of the SHS-PY solution can be avoided by using a 
different version of SHS potential, for which the OZ equations are analytically 
solvable under a different closure, i.e. the "mean spherical approximation" 
(MSA) [8-10]. In this alternative SHS model, the adhesive term is deflned 
starting from a Yukawa attractive tail and taking the "sticky limit" of infi- 
nite amplitude and vanishing range, with their product remaining constant. In 
Ref. I it was pointed out how, in the particular case of factorizable stickiness 
parameters in the adhesive term, the SHS-MSA solution leads to closed ana- 
lytical expressions for scattering intensity and other "global" structure factors, 
which are valid for arbitrary p and therefore are suitable to fit experimental 
scattering data from polydisperse systems. 

In the present work we shall review the main points of Ref. I, adding some 
new results and remarks. We shall then apply the SHS-MSA solution to small- 
angle-neutron-scattering data from a colloidal fiuid formed by sterically sta- 
bilized silica particles. 



2 Model and scattering functions 

Our SHS potential is defined by 

Pu,j{r) = Pufir) - K,^a-^ 5+(r - a,,), (1) 

with factorizable stickiness parameters [11,12] 

Ki^^Y.Y^, >0, (i,j = l,...,p), (2) 
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which mcasTirc the strength of surface attraction. In Eq. (1), u^^{r) is the 
HS potential and the next term stems from the "sticky hmit" of a Yukawa 
tail [13]; /3 = (/cgT)"^, aij = {ai + Uj)/2, with cTj being the HS diameter 
of species i, and S^{x) is the asymmetrical Dirac delta function defined by: 
dx f{x)5+{x — c) — f (c"*") , if a < c < 6. The MSA solution for this model 
is represented by the factor correlation functions 



0, 



elsewhere 



Lij ^ T < CFij 



(3) 



where L^j — {ai — aj)/2 and the expressions of Oj and bi can be found in Ref. 
I. Once that the qij{r) are known, all correlation and scattering functions are 
in principle computable. The partial structure factors are given by [7,15] 

S^Ak)=j:Q;r^.{k)QJ^{-k), (4) 

m=l 



where Qij (k) = Sij — 27r(pjPj)^/^gij (k) , with qij (k) being the unidimensional 
Fourier transform of %(r) {5ij the Kronecker S, pi the number density of 
species i). The scattering intensity as well as other ^'/o&a/ structure factors can 
then be calculated, all these scattering functions being weighted sums of Sij{k), 
like J2ij=i Wi {k) Wj {k) Sij {k) (here, the asterisk means complex conjugation). 
As an example, we mention the measurable average structure factor 

Suik) = E {xiXjf'Fi{k)F*{k)Sij{k) /Y,Xi\Fi{k)\\ (5) 

i,j=l i=l 



where Xi is the molar fraction of species i and Fi{k) its form factor. Eq. 
(4) shows that the main difficulty in the computation of Sij{k) lies in the 
inversion of a p x p matrix Cl{k), which usually becomes a formidable task 
with increasing p. Nevertheless, as discussed in Ref. I, when K^j is factorizable 
as in Eq. (2), Q{k) becomes a dyadic (or Jacobi) matrix, i.e. Qij{k) — Sij + 
Y,^=i ci'i^\k)b^^\k) , which admits analytic inverse for any p (in general, the 
same property does not hold for Baxter's SHS-PY solution [14]). Another 
important consequence of the dyadic structure is that it allows to compute the 
global scattering factors bypassing a preliminary calculation of the individual 
Sij{k). Indeed, when Q(A;) is dyadic the sums Y!i^j=iWi {k)w* (k) Sij (k) can 
be worked out analytically (see Refs. [7,15,16] for details), unlike the most 
common case where they are performed numerically by evaluating p{p + l)/2 
independent contributions Sij{k) at each k. When fitting experimental data, 
the availability of closed analytical expressions directly for scattering intensity 
and Sm (k) clearly represents a great advantage. 
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3 Polydispersity effects in SHS fluids 



3.1 Size and stickiness distributions 

To represent size polydispersity, we select p possible diameters uniformly dis- 
tributed in an interval (0, cXmax) with mesh size Act = 0.02 (a). Each diameter 
(Tj characterizes a different component. Its molar fraction is Xi — f{ai)Aa, 
with / expressed by a Schulz distribution, i.e. 

f{a) = 6V'^-^e-''7r(a) (a > 1), (6) 

where F is the gamma function, (a) the average diameter, a — b = a/ (a) 
and Sa = [(o"^) — (o")^]^^^/ (o") measures the degree of size polydispersity. We 
use p = 85 components when = 0.1, and p = 175 when s^- = 0.3 . 

The model can in principle be worked out for any choice of Yi, which meets 
the requirements that Yi must be lenghts and Kij = YiYj must be proportional 
to /? = {kBT)~^. However, two choices are particularly interesting: 

Model I (polydisperse in size but monodisperse in stickiness), which assumes 
that all particles have the same stickiness, i.e. 

y, = y = 7oM, with ^i^^^l-. (7) 

Here, 70 measures the strength of surface adhesion; £0 is an energy, and T* 
a reduced temperature. High T* (low 70) values correspond to weak surface 
adhesion or high temperature; vice versa, low T* (high 70) values mean strong 
attraction or low temperature. 

Model II (with stickiness polydispersity linearly related to the size one). Since 
it is reasonable to expect that larger particles attract each other more strongly, 
a simple way of introducing stickiness polydispersity is the choice 

Yi ^joai. (8) 



3.2 Regimes for SMik) and generalized Boyle temperature 

Polydispersity effects on Suik) can be best described after recalling some 
features of the monodisperse SHS structure factor. 
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In (Smono (k) the presence of strong adhesive forces may be revealed mainly 
from its behaviour in the low-A; region. Here. 5*,^^^^ (k) may differ significantly 
from 5*^0110 (k) , exhibiting a drastic increase as — > (see Figure 1 for s^. = 0). 
More precisely, in Ref . I we recognized the existence of two different "regimes" 
for (k) respectively above and below the Boyle temperature, — 3, 

where attractive and repulsive forces balance each other in such a way that 
B2 vanishes (the second virial coefficient B2 appears in the low-density expan- 
sion of S^ono (0) = 1 - 2B2{T*)p + 0{p^) ). When T* > {B2 > 0) the fluid 
behaves like pure HS without stickiness, repulsive forces are dominant, and 
'S'mono (0) - which is related to compressibility and density fluctuations - de- 
creases with increasing volume fraction rj. When T* < {B2 < 0) S'mono (0) 
has a non-monotonic dependence on 77. After an initial increase (which may 
become very strong, signalling the approach to a gas-liquid phase transition 
with critical temperature T* ~ 1-61), an inversion occurs and S'mono (0) de- 
creases with rj. In other words, below Tg the attractive forces are dominant at 
low T), whereas repulsion prevails at higher 77. 

For polydisperse SHS, there exist two similar "regimes" for Su {k). The new 
feature is that a crossover between these regimes can now be induced not only 
by a change of T*, but also by a change of polydispersity Sa- (at fixed T*). 

Figure 1 illustrates the polydispersity effects on Suik), displaying its be- 
haviour at T* = 2.04 (70 = 0.7) and rj = 0.2, as increases from (monodis- 
perse limit) to 0.1 and 0.3. We employed form factors Fi{k) for homogeneous 
spherical scattering cores inside the particles, with scattering core diameters 
^scatt coincident with the HS ones, i.e. o-|=^t* = o-,.. Polydispersity progressively 
dumps all oscillations in the first peak region and beyond. Here, Model I and II 
are nearly equivalent, whereas their Suik) may strongly differ near the origin 
as a consequence of the different stickiness distributions. In fact, increasing 
(at the fixed T* value) does not produce marked differences in Model II. On 
the contrary, in Model I a large polydispersity can overwhelm the attractive 
effects, leading to a strong HS-like decrease of S'm(O) and thus preventing the 
possibility of an instability in the system. 

In Ref. I this scenario was interpreted in terms of a generalized Boyle tempera- 
ture, Tg p. From the low-density expansion of Sm(0) = 1-2B2^f(T*)p + 0{p^), 
we derived i?2,F as a generalization of B2 when all form factors are included. 
Then B2^f{T^p) = O defines T^p{sa-), which is a decreasing function of s^- 
(its expression is given in Ref. I). Clearly, T^^p{sa = 0) = (TB)mono = 3. In a 
(so-, T*) diagram, above the curve T* — T^^pisa) there is the region correspond- 
ing to a "HS-like, repulsive regime" ; below, there is the "attractive-regime" re- 
gion. Now, as shown in the inset of Figure 1, Tg p{s„) of Model II decreases very 
slowly, asymptotically approaching 2.57; therefore, all three (so-,T* = 2.04)- 
states of Figure 1 lie in the lower region. On the other hand, T^^p^s^) of 
Model I is a rapidly decreasing function, with T^^p{sa — 0.1) = 2. '79 > T* 
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and Tj^ p(.Scr = 0.3) = 1.68 < T* . This argument interprets the crossover from 
"attractive" to "repulsive" regime in terms of an increasing polydispersity. 



4 Fit of experimental data 

To investigate polydispersity effects in real colloidal fluids, we now apply our 
Models to experimental data for silica particles coated with octadccyl chains 
and dispersed in benzene. For this system Duits et al. [3] published a wide 
set of small-angle neutron-scattering results, which exhibit a crossover from 
"HS-likc repulsive" to "attractive" regime controlled by variation of T: on 
lowering the temperature, attractive effects increase and a reversible phase 
separation into two phases of different density occurs at t = 32.5 °C. Duits 
et al. [3] fitted their experimental data using an implementation by Robertus 
et al. [1] of Baxter's SHS-PY model, with stickiness independent of size and 
size polydispersity represented by a rather limited number of components (up 
to p = 9), due to the aforesaid shortcomings of the PY solution. The form 
factors Fi{k) were calculated with a model of spherical three-layer particles, 
and a difference between HS and scattering core diameters, with cTj > crf^**", 
was admitted. 

Since our aim is not an exhaustive analysis of all data by Duits et al. [3], 
in Figure 2 we consider only two representative cases, with r] = 0.28 and 
t — 51.6 °C and 35.4 °C, respectively. The Suik) of the first sample (Figure 
2a) is typically "HS-hke", whereas at low k the curve of Figure 2b has an 
upswept shape characteristic of strong attraction. For simplicity, we calculate 
the form factors with a model of homogeneous spherical scattering cores. We 
have checked that this choice is essentially equivalent to the three-layer one, 
if appropriate average contrast and diameter ( (cr^'^^**) = 39.8 nm ) are used, 
as suggested by Ref. [3]. 

In Figure 2a we first fit the experimental data by employing the monodis- 
perse SHS-MSA model (all particles having the same size and stickiness; 
So- = 0), with two free parameters, T* and a > (cr^'^^**) . The best fit results 
are: (T*, a) = (4.2, 45.6), with the value of T* confirming that the system 
is in the "repulsive" regime. Clearly, this model overestimates the first peak 
height and does not reproduce the rising part of the curve adequately. Both 
these drawbacks can be overcome by taking size polydispersity into account. 
Indeed, both Model I and II yield excellent agreement, and their results are 
nearly identical, indicating that the choice of stickiness distribution is irrele- 
vant in this HS-like case. To use only two free parameters, i.e. (a) > (o'^'^^**) 
and Scr, for Model 1 and 11 we fix T* on the value found with the monodis- 
perse model. The best fit results are: ((cr) , s^) = (44.9, 0.12) for Model I, and 
(44.2, 0.15) for Model II (using three free parameters, i.e. T*, (cr) , So-, does 
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not alter oTir qualitative results). We remark that both our estimates of Sa- are 
consistent with the value 0.11 of Ref. [3]. Moreover, comparison with Figure 
6 of Ref. [3] shows that Model I improves over its PY counterpart of Duits et 
al. This may be due to our use of a much larger number p of components to 
represent the size distribution. 

In Figure 2b we use a similar strategy. Here, however, we have found impossible 
to satisfactorily fit the entire experimental curve with the monodisperse model. 
Hence, we have forced T* to the value 2.3 which best reproduces the data 
within the restricted /c-range (0.03, 0.15) nm"-*^ (note that this value of T* 
indicates that the system is in the "attractive" regime). Then, upon fitting 
the complete set of data with one free parameter, we find (a) = 43.0 nm for 
the monodisperse model. On fixing T* = 2.3, the best results for Model I and H 
are: {{a) , So-) = (41.9, 0.12) and (41.5, 0.13), respectively. Once again we note 
the importance of the polydispersity effects, which lower the first peak height 
correctly and improve the general agreement. Both Models are satisfactory, 
although Model II appears to be slightly superior, indicating, perhaps, that the 
difference of stickiness distribution is already becoming significant. However, 
this conclusion cannot be definitive. In fact, further calculations (not reported 
here) under different conditions, yield ambiguous results. Moreover, Figure 1 
shows that the differences between Model I and II are small at low values, 
as unfortunately occurs in the colloidal system investigated by Duits et al. [3]. 
On the other hand, our results suggest the opportunity of further experimental 
investigations on similar fluids with larger polydispersity, where an analysis of 
Suik) at low k could discriminate between different SHS Models and provide 
useful information about the unknown real stickiness distribution. 



5 Conclusions 

The present study on a new polydisperse SHS model has once again shown 
the importance of taking polydispersity effects into account when experimen- 
tal scattering data are analyzed. We have found that even a small degree 
of size polydispersity, — 0.1, can strongly modify Suik) with respect to 
the monodisperse case. We have pointed out the existence of "repulsive" and 
"attractive" regimes for Su{k), with a crossover interpreted in terms of a "gen- 
eralized Boyle temperature" , dependent on the degree of polydispersity. In the 
presence of polydispersity, a crossover between these regimes can induced not 
only by a change of T*, but also by a change of (at flxed T*). 

In addition to size polydispersity, we have considered two simple stickiness 
distributions (Model I and II). Our results indicate that the particular choice of 
stickiness distribution is irrelevant in fluids with small size polydispersity and 
weak coupling (i.e. weak attraction and high temperature). On the contrary. 
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stickiness polydispersity effects can be revealed from tlie \ow-k beliaviour of 
Suik) when size polydispersity is large and coupling is strong (strong adhesive 
forces and low temperature). 

As an application, we have employed Model I and II to analyze experimental 
scattering data for a colloidal fluid with short-range attractions. Our best fit 
results are satisfactory and suggest the importance of the stickiness polydisper- 
sity (Model II). However, the considered experimental system lies in a regime 
where it is impossible to unambigously discriminate between our two Models. 
New scattering studies on fluids with larger size polydispersity and stronger 
adhesive forces would be desirable to get more insight into this problem. 
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Fig. 1. Structure factor Suik) of the SHS-MSA Models I and II at T* = 2.04 and 
r) = 0.2, for several degrees of polydispersity (the curves of Model I are shifted 
upwards by 1.5 units to avoid overlapping). In the inset, the solid curves represent 
Tg p(scr) of the two Models, the dashed line corresponds to the isotherm T* = 2.04 
and the black squares on it are the considered states. 
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Fig. 2. Fit of experimental neutron scattering data for Suik) at ij = 0.28 and two 
different temperatures (taken from Ref. [3]), by using monodisperse and polydisperse 
SHS-MSA models (the curves of part (a) are shifted upwards by 2.5 units). 



9 



